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Abstract. In line with the stability theory of continuous dynamical 
systems, Lyapunov exponents of cellular automata (GAs) have been con¬ 
ceived two decades ago to quantify to what extent their dynamics changes 
following a perturbation of their initial configuration. More precisely, 
Lyapunov exponents of GAs have either been understood as the rate by 
which the resulting defect cone widens as these dynamical systems are 
evolved from a perturbed initial configuration, or as the rate by which 
defects accumulate during their evolution. The former viewpoint yields 
insight into the extent of the affected region, whereas the latter tells us 
something about the intensity of the defect propagation. In this paper, 
we will show how these viewpoints can be united by relying on Lyapunov 
profiles of GAs. 


1 Introduction 

Ever since their postulation by von Neumann [9], researchers have often been 
struck by the intriguing spatio-temporal dynamics that cellular automata (CAs) 
can bring forth in spite of their intrinsically simple nature [4,10]. In order to 
arrive at a comprehensive analysis of their behavior, different measures have 
been proposed, such as Lyapunov exponents [2,6] and others [5]. With respect 
to Lyapunov exponents, two distinct viewpoints have emerged during the last 
two decades. The first one traces back to a suggestion by Wolfram [10], which 
was formalized by Shereshevsky [6], and later on adopted by several authors [3, 
7], and which identifies Lyapunov exponents of elementary CAs with the rates 
by which the damage front moves to the right and to the left. Consequently, 
every elementary CA is characterized by means of two exponents, which obvi¬ 
ously might differ. On the other hand, the viewpoint put forward by Bagnoli 
et al. [2] involves only one Lyapunov exponent for every CA that quantifies the 
exponential rate by which defects accumulate if the CA is evolved from a initial 
configuration with a single defect without accounting for the direction in which 
defects propagate [2]. 
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The definition of Lyapunov exponents by Bagnoli et al. [2] actually consti¬ 
tutes the counterpart of the one by Shereshevsky [6] because accumulation is 
neglected while directionality is considered in the latter, whereas the opposite 
holds for Lyapunov exponents according to Bagnoli at al. In other words, the for¬ 
mer viewpoint gives insight into the intensity of the defect propagation, whereas 
the latter indicates how fast an initial perturbation spreads as the CA evolves 
and therefore how fast the defect cone widens as the CA evolves. Clearly, each 
viewpoint has its advantages, but uniting them in a single concept enables a 
comprehensive stability analysis of CAs. 

In this paper, we will show how so-called Lyapunov profiles of CAs constitute 
such a uniting framework that allows for assessing both the intensity of the 
damage propagation and the growth rate of the defect cone. 

2 Lyapunov profiles of cellular automata 

2.1 Preliminaries 

In the framework of this paper, we will restrict ourselves to elementary CAs 
which can be conveniently represented by a quintuple = {T, S, s, N, f2). 
Therein, T denotes a countably infinite tessellation of a 1-dimensional Euclidean 
space, consisting of consecutive intervals Ci, i G N, typically referred to as cells, 
and S' is a set of two states, here S = {0,1}. Further, the output function 
s : T X N —>■ S gives the state value of cell Ci at the t-th discrete time step, N{ci) 
is the ordered list of neighbors of Ci, i.e. N(ci) = (ci_i, Ci, Ci+i), and finally the 
transition function 17 : N —)■ S that governs the dynamics of each cell Ci reads 

s{c^,t+ 1) = n s(c^,t), s{ci+i,t)) . 


2.2 Rationale 

The Lyapunov exponent of a CA according to Bagnoli et al. [2] is computed by 
tracking all defects that emerge from the introduction of a single defect in the 
CA’s initial configuration, i.e. by flipping one of the cells’ state, as follows: 

A = lim ^log V (1) 

t-J-oo t \£0 / 

where et denotes the total number of defects at the t-th step and a defect is 
defined as the smallest possible perturbation of the CA, such that it may be 
conceived as an object that flips that state of the cell in which it resides. Upon 
introducing such a defect to one of the CA’s cells, it will propagate and affect 
increasingly more cells because the dynamics of the CA is governed by its neigh¬ 
borhood configurations. The quantity et should be computed by evolving both 
the initial configuration sq and its perturbed counterpart Sg, for which it holds 
that d(so,s5) = I, for one time step, after which a replica s*(-,l) of s(-,I) is 
created for every cell Ci for which it holds that s(ci, I) ^ s*(ci, 1) in such a way 
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that d{s{-, 1), s* (•, 1)) = 1. Then, all replicas are evolved one more time step and 
the resulting configurations s*(-,2) are again compared with s(-,2), such that a 
new set of replicas can be constructed, after which this procedure is repeated 
until the CA evolution halts. So, when computing the total number of defects 
up to a given time step, both the position and multiplicity of defects is known, 
but this information is lost when computing the Lyapunov according to Eq. (1). 

Clearly, if we redefine e* as a vector whose f-th element ej is the number 
of defects in cell Ci at a given step, we have a construct that incorporates both 
the multiplicity and position of defects. For comprehensiveness, it should be 
mentioned that Bagnoli et al. [2] do incorporate such a damage vector to define 
the Lyapunov exponent of elementary CAs, but they discard its information 
content as their further analysis is based upon the sum of its elements. Besides, 
they call the quantity given by Eq. (1) the maximum Lyapunov exponent (MLE) 
of a CA, by analogy to continuous n-dimensional dynamical systems where the 
terms refers to the largest exponent in the Lyapunov spectrum that contains the 
perturbations, and rates of separation, in all n different directions. By contrast, 
in the case of binary CAs that Bagnoli et al. [2] consider, they define only one 
exponent, thus the meaning of the term ‘maximum’ in their definition is unclear. 

In this paper we interpret a CA on an array of n cells as an n-dimensional 
system, thus it may be anticipated that a CA evolution generates a spectrum 
of Lyapunov exponents, one per each cell of the array. The largest of these is 
naturally called the MLE. As et quantifies how the number of defects grows in 
each possible direction, i.e. in each cell, and therefore provides insight into the 
damage propagation resulting from a single defect in each direction, we proceed 
to take a closer look at what exactly we might learn from . 


2.3 Lyapunov profiles 

At this point, we shall define the finite-time Lyapunov profile P of a CA as 

At = ^ log (ct) 

where log is applied element wise. Hence, the Lyapunov profile consists of the 
time-averaged expansion rates in each of the possible directions. The elements 
of At may be understood intuitively as the time-averaged exponential rates by 
which the number of defects grows in the cells of the CA, and we could look at the 
largest among them as the MLE of the CA. Besides, the width of the defect cone 
is given by |{i | loge^ ^ this point it should also be mentioned that 

the defect distribution does not translate directly to the distribution of damages 
in the CA configuration space. In the remainder, we will refer to ^ log(ey) as 
the f-th finite-time Lyapunov exponent of a CA, denoted A^. 

Given the fact that the propagation of defects is constrained by the size of 
the neighborhood N{ci) [2], there must exist an upper bound on At’s elements. 
Recalling that in the worst-case scenario a defect in a cell Ci will propagate to 
each of its three neighbors if the CA is evolved for one time step, and this occurs 
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for all cells in T, it is easy to see that the upper bound on At is given by the 
T + 1-th row of the trinomial triangle in such a way that 



where fg denotes the index of the cell that was perturbed at t = 0. The sum of 
its elements on the T + 1-th row are 3^, which is the upper bound on the total 
number of defects at T, which thereby indicates why the MLE according to [2] 
can be at most log(3). 

The Lyapunov profiles of rules 105 and 150 at T are completely defined by 
the T -I- 1-th row of the trinomial triangle because all cells in these CAs pass on 
defects to all their neighboring cells. Besides, their maximum will be given by 


maxAr = 




or equivalently [1]: 


max At 



( LWJ i t-2o 

i+ui: If 


( 2 ) 


y i=i fe=i j 

which provides us with a means to normalize the Lyapunov spectra of elementary 
CAs, such that a mutual comparison becomes possible. 


3 Results 

3.1 Experimental setup 

In the remainder of this paper we will focus our attention on the 88 minimal 
elementary CAs as defined in [8]. For each of them, the propagation of de¬ 
fects emerging from a single defect was tracked for 5000 time steps in a one¬ 
dimensional system consisting of 10001 cells, as such mimicking a system of 
infinite size. In order to enable a meaningful comparison between the Lyapunov 
profiles of the considered CAs, their elements were normalized by means of the 
upper bound on the exponential propagation rate as given by Eq. (2). Although 
the Lyapunov profiles will depend in in many cases on the initial configuration 
from which they are evolved (Class 2 and 4 rules), we restrict ourselves in this 
preliminary work to uniformly chosen random initial states and a single initial 
defect. A way to eliminate some of the dependence is to replace the singleton 
with a large zone of defects, which makes it very unlikely that defects die by 
chance. 
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3.2 Lyapunov profiles of elementary cellular automata 

The benchmark rule As a kind of benchmark, Fig. 1 shows the normalized 
Lyapunov profile of rule 150, which is equivalent to the one of rule 105. For 
comprehensiveness, the cell indices along the horizontal axis are positioned in 
such a way that index 0 refers to the cell where initial defect was introduced. 
As anticipated when deriving the upper bound on the MLE, the normalized 
profile of this rule reaches a maximum 1 in the initially perturbed cell. The 
discrepancy between the line Ay = 1 and the profile gives an indication of the 
delay in the defect propagation that is imposed by the finite speed of information 
transmission in CAs. Indeed, it takes t time steps before a cell that is t cells 
apart from the initially perturbed one can start propagating defects. As such, if 
one would normalize Ay’s elements with respect to the corresponding trinomial 
coefficient, and not with central trinomial coefficient of the T + 1-th row of 
the trinomial triangle, such a normalized Lyapunov profile would result in a 
horizontal Ay = 1. At this point, we opted not to normalize the Lyapunov 
profile in this way because one would then somehow lose the information on the 
delay of the defect propagation that is imposed by the finite speed of information 
transmission in CAs. 



Fig. 1. Normalized Lyapunov profile of the ECA rule 150 that was obtained after 5000 
time steps for a system of 10001 cells. 


Class 3 rules Figure 2 depicts the normalized Lyapunov profiles of some of the 
most well-studied ECA rules that belong to Wolfram’s class 3 [10]. It is clear 
that the profiles of ECA rules 30 and 90 lie below the one of rule 150, but there 
are pronounced differences between the profiles of these Class 3 rules. The one of 
rule 90 shows the most agreement with the one of rule 150, with a very similar 
shape, but with a maximum of approximately 0.6 instead of 1. It is importance 
to notice that it has positive values for all cells, which means that the defect 
cone widens at maximum speed, just as in the case of rule 150, but its intensity 
of defect propagation is lower. The defect cone for rule 30 is much smaller as it 
exhibits a sharp transition at i « —2300 from expansive cells where Ay > 0 to 
cells that are not yet affected by defects, and hence have Ay = —oo. Moreover, 
the one of rule 30 is asymmetric, in the sense that its damage cone propagates to 
the right at maximum speed, whereas it travels to the left at a speed lower than 
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half of this maximum and its maximum is not located at the initially perturbed 
cell. 



(a) Rule 30 (b) Rule 90 


Fig. 2. Normalized Lyapunov profiles of some Class 3 EGAs that were obtained after 
5000 time steps for a system of 10001 cells. 


Class 4 rules Figure 3 depicts the normalized Lyapunov profiles of two Wol¬ 
fram’s [10] Class 4 EGA rules that are included among the 88 minimal EGAs. 
Compared to the Lyapunov profiles of the exemplary Class 3 rules, the ones of 
the Class 4 rules are more intriguing in the sense that they typically exhibit 
one or two phase transitions, i.e. sharp transitions from expansive cells where 
> 0 to not yet affected cells (rules 106 and 110), are asymmetric (rules 106 
and 110), and/or non smooth (rule 106). Besides, just as we observed in the case 
of rule 30, the maximum of the Lyapunov profile does not necessarily occur at 
the cell where the defect was introduced at the beginning of the evolution. This 
skewness is especially pronounced in the case of rule 110. 



(a) Rule 106 (b) Rule 110 


Fig. 3. Normalized Lyapunov profiles of the Class 4 minimal EGAs that were obtained 
after 5000 time steps for a system of 10001 cells. 
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Class 2 rules Although most non-trivial Lyapunov proHles are found for rules 
belonging to Wolfram’s Classes 3 and 4, it should be emphasized that also Class 2 
rules can give rise to non-trivial Lyapunov profiles. More precisely, such rules 
either give rise to a highly discontinuous profile that are — oo, except in a few 
cells, or a truly meaningful profile. For instance, within the family of minimal 
ECAs, rules 28, 33, 37, 73, 108 and 156 give rise to such discontinuous profiles, 
which indicate that defects can accumulate as the ECAs evolve, but they remain 
localized and cannot propagate beyond certain limits. We opt not to include 
exemplary profiles of these rules because they are — oo almost everywhere. On 
the other hand, there are a few rules, namely 6, 57 and 62, which give rise to 
smooth prohles, similar to the ones of the Class 3 rules (Fig. 4). Other class 
2 rules within the studied ECA family lead to a Lyapunov profile that is —oo 
everywhere. 

As indicated by the profile of rule 57, even Class 2 rules can give rise to 
a Lyapunov profile that stretches across the entire light cone, i.e. [—T, T], but 
this is the only Class 2 rule that exhibits such a pronounced defect propagation. 
For the other Class 2 rules that give rise to a non-trivial Lyapunov profile, it 
is typically lower and narrower than the ones of Class 3 and Class 4 rules. The 
peculiar Lyapunov profile of rule 62 can be better understood by looking at its 
corresponding heat map, which displays the time evolution of the normalized 
Lyapunov profile (Fig. 5). This clearly shows that the defect propagation is 
confined to a region of which the right side extends at maximum speed, whereas 
its left side tends to be pushed inwards substantially from time to time. 



(a) Rule 6 (b) Rule 25 (c) Rule 38 


Fig. 4. Normalized Lyapunov profiles of ECA rules 6, 57 and 62 that were obtained 
after 5000 time steps for a system of 10001 cells. 


Conclusions and future work 

In this paper it was shown how the distinct viewpoints on Lyapunov exponents 
of CAs can be united by considering their so-called Lyapunov profiles. These 
prohles consist of the time-averaged rate of defect accumulation in each of the 
possible directions that is allowed by the CA’s structure. As a next step, we will 
focus on whether rigorous results can be obtained for the Lyapunov prohles and 
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Fig. 5. Heatmap of Rule 62 that was obtained after 500 time steps for a system of 1001 
cells. 
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the therewith defined MLE, their extension to arbitrary two-state families, as 
well as their formulation in the case of multi-state CAs. 
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